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ABSTRACT 

We consider the one-dimensional Hubbard model with the infinitely strong 
repulsion. The two-point dynamical temperature correlation functions are cal- 
culated. They are represented as Fredholm determinants of linear integrable 
integral operators. 
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1 Introduction 



The one-dimcnsional Hubbard model ^ is one of the most interesting and im- 
portant model of strongly correlated electrons (see, e.g., book [^). The hamilto- 
nian of the model with the chemical potential h in the external constant field B 
describes interacting fermions on the one-dimensional periodic lattice of length 
L: 

L 

L L L 

-h^irix.,-] + n^,i) +B^{nx,^ - n^.i) + U^n^.-^n^^i. (1) 



x = l 



Here „ and 4^ (x = 1, . . . , L; a =T, U Cq,, 
fermion fields with the anticommutation relations 



are canonical 



(2) 



The density operators of electrons with spin up (a =t) and down {a =i) are 
'^a; Q = 4 a- The solution of this model by means of the two-component 
(nested) Bethe Anzatz ^ was given in 

The temperature correlation functions are defined in a standard way as tem- 
perature normalized mean values. In the present paper we consider two-point 
correlation functions of canonical field operators. 
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where a, (3 =t, [. The dependence on time is introduced as usual. 
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The trace Sp is taken in the Fock space T of dimension 4-^ where the hamihonial 
of the Hubbard model acts. The correlation functions are mostly interesting in 
the thermodynamic limit (i — > oo; h and B are kept fixed): 

(...)m= lim(...)(^'^). (6) 

L — ^oo 

Our aim is to present the results on calculation of the correlation functions 
in the Hubbard model in the strong coupling limit, U — > +00. The two-point 
temperature correlation functions arc represented as Fredholm determinants of 
linear integral operators of a special kind (the "integrable" integral operators, 
in the sence of the paper |^). 

It is to be mentioned that the recent progress in calculating correlation func- 
tions of quantum solvable models is based on the fact that they are governed 
by classical integrable differential equations. That the language of classical dif- 
ferential equations is quite natural for the description of quantum correlation 
functions was realized a time ago ||, ||, |l0j . The idea of the approach suggested 
in [0 H is to consider the Fredholm determinant in the representation for a 
correlation function of a quantum integrable model as a tau-function for a clas- 
sical integrable system (see also the book where the results for the simplest 
model of one-dimensional impenetrable bosons are reviewed). The necessary 
first step (which is also of interest by itself) in this approach is to represent 
the correlation function as the Fredholm determinant of a linear "integrable" 
integral operator. The first determinant representation of this kind was given 
in ^1 for the equal-time temperature correlators of the onc-dimcntional im- 
penetrable bosons. For the two-component one-dimensional impenetrable Bose 
and Fermi gases the representations of this type were obtained recently in papers 

The Hubbard model is a lattice analogue of the two-component Fermi gas. 
It should be noted that the physics of the Hubbard model at [/ = -l-oo has some 
interesting properties (see, e.g., ||l^). In particular, the ground state becomes 
degenerate at zero external field, i.e., the point S = is the point of the phase 
transition (the same take place in the Fermi gas as well). 

We use the technique of calculating temperature correlation functions of 
exactly solvable quantum two-component models with infinitely strong coupling 
developed in papers 17 for one-dimensional two-component Bose and Fermi 



gases. The starting point of the approach is using the eigenfunctions of the XX^ 
spin chain with the periodic boundary conditions to solve the auxiliarly lattice 
problem of the nested Bethe Ansatz. This gives an explicit expression for the 
eigenstates of the two-component models at the point of infinite coupling. 



2 The hamiltonian and eigenstates at t/ = oo 

In the strong coupling limit, the states with double occupancies (i.e., the states 
with a pair of electrons at least at one cite of the lattice) have infinite energy 
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and therefore are absent in the physical space of states. The dynamics of the 
Hubbard model in this limit is described by an effective hamiltonian which can 
be written in the form 119, M 



H = V 



where 



x = la = ],i 
L L 

X—1 X—1 



V, 



(7) 



p=n(i 



)• 
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The projector V extracts the physical space of states Ti of dimcntion 3^ from 
the space J- (of dimention 4^) in which the canonical Hubbard operators act. 

Eigenstates of the effective hamiltonian H (all belonging to the space TL) 
can be constructed analogously to the case of the two-component impenetrable 
Fermi gas |l^ . The eigenstates have the form [|| 



zi,...,zn — 1 ai,...,aiv — f,J, 
cL,a,---4„,a„|0), (9) 

where the Fock vacuum |0) is defined as usual, c, JO) = 0, (0|4 a = 0, (0|0} = 1. 
The eigenstates are parametrized by set of A'' unequal real numbers (quasimo- 
menta of N particles), k = fci, . . . , fcjv, and the set of M unequal real numbers 
(quasimomenta of the auxiliarly lattice problem), A = Ai, . . . , Am- The num- 
ber A'' is the total number of electrons in the state while the number M is the 
number of electrons with spin down. Thus M of values a in (||) are J,, and 
N — M are f. The momenta k and momenta A are not arbitrary but must obey 
the system of the Bethe equations (the nested Bethe Ansatz) H, ^ ||]. In the 
case U = +0O using the XXO basis for the auxiliarly problem results in the 
simplified form of the Bethe equations. 



a — 1, . . 

b- 



1,. 



where 



M 
6=1 



(10) 
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This system can be solved explicitly; the equations are separated by the substi- 
tution 



A 
L' 



(12) 



The permitted values of the momenta are 

T 

2n 

iV V 2 



{ka)j 



-J, j ^0,...,L-1, 

N 1 + {-l)N-M 



+ 1 



I = 0, 



,N-1. 



(13) 



The wave function xn,m is 
X%\ir"{zi,...,ZN\k;X) 



1 

m 



detAr{e*'''"^^} 



(14) 



The sum here is taken over the permutations of N numbers, P : (1, . 
(Pi, . . . , P/v). The function 9{zi < . . . < z^) is equal to 1 if zi < . . . < za? and is 
equal to zero otherwise. The determinant detArje*'^"^''} denotes the determinant 
of the N X N matrix with elements e''^"^''. The spin part of the wave function 
is described by components of the 2^-dimensional vector \£,n.m{X))- As in the 
case of the two-component impenetrable gas, we use the vectors \S,n.m{X)) to 
be eigenvectors of the hamiltonian of the XXO spin chain in the form given in 



N 



(15) 



ni,...,njif = 1 

with the wave function 



1 

M! 



^NMi"^!^ ■ ■ ■ :"Af |A) 

Pauli matrices are defined as usual 



l<j<l<N 











r 









detM{e'^°"^}. (16) 



(17) 



The vacuum | f|-jv) of the spin chain with N sites is defined as cr^^| f^Ar) — 0, 
(^^1%) =0, (^atI^jv) = 1, i.e., 

IM = ®:Lifn') • (18) 
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Eigenvalues of the hamiltonian (|^) on the eigenstates (Q), 

H |*Ar,Af(fc, A)) = EN^Mik) \^n,m{K A)), (19) 

are 



N 



EN,M{k) = Y.^{ka) -hN + BiN - 2M), (20) 



a=l 

where e{k) is one-particle dispersion 

e{ka) = -2 cos ka = -2 cos ^fc^ + . (21) 

Thus, the eigenenergies depend on Ai,...,Am via their sum, A. It is used 
essensially in our approach. 

The states \'^N.M{k; A)) form a complete orthogonal set in the physical space 
Ti. (of dimention 3^); the normalization is 

(^AT^Hfc; X)\'fN,M{k; A)) = L^'N^'. (22) 

The effect of the infinitely strong interaction consists in contraction of ob- 
servables (quantum operators) onto the space Ti.. The unit operator in Ti is 
just the projector V in T. Due to the completness, this operator possesses the 
decomposition 

r^Y. l*^.M(fc;A)}-^^^(*jv,M(fc;A)|, (23) 

states 

where the summation is performed over all 3^ basis states in Ti., 

L N 

E-EE E • (24) 

states Af=0 Af=0 '=l< - <fcjv 
Xi<...<Xm 

Due to (psf), the effective hamiltonian H can be also defined by its spectral 
decomposition 

H^J2 \^N,M{k;\))^gr0-('^^Mk;\)\. (25) 

states 

The temperature correlation functions in the strong coupling limit acquire 
the form: 



(4,a(i)Co.a(0)) 
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(c.aWcL(o)) 
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Here the trace Tr (which is different from Sp in (H) and 
the space Ti., i.e., by definition, 



Tr 



E 



1 



(*7v,Af(fc;A)|---|^'7v,Af(fc;A)}. 



(26) 



, is to be taken in 



(27) 



The tildes over the Hubbard operators in the numerators in (^ 
contraction of these operators onto the space 7i, i.e.. 



V 



the 



(28) 



The time dependence is described by the relations similar to (|^) with the ef- 
fective hamiltonian H given by (Q). The equations (^6|) exhibit the recipe for 
calculation of the correlators in the strong coupling limit. Indeed, all states be- 
longing to the subspace T\H (of dimention 4^ — S'^) in the decomposition of the 
unit operator in the Fock space T inserted between canonical field operators in 
the numerators in (^) do not contribute to the correlators due to infinitely strong 
oscilation for any small (but finite) value oi t as U tends to infinity. The same 
happens with the correlators of density operators, but since [P , nx^a{t)\ = such 
intermediate states do not contribute already due to the temperature exponen- 
tial. Thus, equations ( p6|) together with ( p8| ) and (|2^) allow one to express the 
correlators in the strong coupling limit through matrix elements (form factors) 
of Hubbard operators ^ [t) , cj, (t) , „ (0 between two states from TL only. 



It is to be emphasized that equations (|2^) will produce the correct answers also 
for the equal-time correlation functions at [/ = oo (since t should be put equal 
to zero after taking the limit U oo). 

The equal-time correlators of the canonical field operators (given by the 
first and the second equations in (^) at t = 0) should satisfy simple relations 
governed by the anticommutation relations between ^ and p. lix ^ y then 

Cj, „ and ^ anticommute, while at x = y and a — (3 the anticommutators are 
different from ^ being equal to 

V 



V 



In particular, it follows from these relations that 



(and the same for , cj^ , ). Having in mind that cj, j^c^. ^ = 



(29) 

(30) 

and that 



V play the role of the unit operator in the space Ti, it is easy to see what the 
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relations between correlators at i = 0, a; = are. Let us note that these relations 
are different from those which one gets from (^) and (H) . It means that the limits 
/7 — > oo and t — > do not "commute" . 

The operators c^_^ and cj. due to (p3|), can be also defined by means of 
their decompositions, e.g., one has 



L-l N 

EE E E 

N=OM=a «i< - -Kgjv fci<. .<fcN + i 

Hl<.-</JAf Ai<...<Aj,,, 



' L^^+^N'^^N + iy 



where (g; /i|fc; A) arc matrix elements (form factors) of the operator 

I between two states from H: 

^l^+'^''''(g;//|fc;A) = (vl/^,Af(9;M)|c,,Tl«'iv+i,Af(fc; A)). (32) 
Explicit calculation results in the following representation 

^iT^''^ /^Ifc; ^) = detM det^+ii? 

{Af+l W ~j 

ix ka - ix'^qb > , (33) 
a=l fc=l J 

where A = J^aLi and 9 = J^tLi t^-b- Elements of the {N + I) x {N + 1) 
matrix are (a = 1, . . . , iV + 1; & = 1, . . . , iV) 

{DU = cot(^^) , {DUr,+i - 1, (34) 

and elements of the M x M matrix W are {a,b — 1, . . . , M) 

N 

(W^),fc = ^e™(^»-'^''). (35) 



n=l 



The similar representation can be obtained for the form factors of operator ^. 
The form factors of operators cj. ^ can be obtained from the form factors of 
operators ^ by means of complex conjugation. The form factors of density 
operators n^.a can be also calculated; the corresponding results are given in our 
recent paper |^ . 

Correlation functions ( ^ ) can be expressed, due to (|3l|), as the sums over 
all intermediate states of squared modula of the corresponding form factors. 
The resulting expressions for the correlators are quite similar to those obtained 
earlier for the two-component impenetrable gas [l^ . On a finite lattice, they 
are rather bulky, being considerably simplified in the thermodynamic limit. 
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3 Correlation functions in the thermodynamic 
Umit 



In this Section the results of our calculation of the two-point temperature cor- 
relation functions in the Hubbard model at C/ = oo in the thermodynamic 
limit, L oo, are given. To derive the representations below from the repre- 
sentations on a finite lattice we use the technique described in detail (for the 
two-component gas) in paper [^ . In the limit, the sums over the intermediate 
states in the expressions for the correlators on a finite lattice can be reduced to 
expansions of the Fredholm determinants of some linear integral operators. 

Thus, the correlation functions considered are represented as Fredholm de- 
terminants of linear integral operators. The integral operators which enter the 
representations for the correlators act on arbitrary function /(fc) according to 
the rule: 

TT 

{A- f){k)= I dk' A{k,k')f{k'), (36) 



where the function A{k, k') is the kernel of the integral operator A. 

In order to write down the representations for the correlators, let us introduce 
some notations. Define the functions 

TT 

r ^—ite{q)+ixq ^—ite{k)+ixk 



E{k)^E{k;x,t) J dq 



tan(2^) 



G = G{x,t) = / dge-'*"(«)+*"«. (37) 



The distance x is an arbitrary integer; the dispersion is e{q) — —2 cos q (see (pl[)) 
where the momenta q and k in ( ^7| ) take any values on the interval [— tt, tt]. In 
particular, G{x,t) = i^Jx{2t) where is the Bessel function {G{x,0) = Sx,o)- 
The following pair of functions play an important role; 

e-{k) ^ E^{k)^/d{k), (38) 
where rj G [— 7r,7r] is a parameter and i?(fc) is the Fermi weight, 

g-B/T 

= 2coshf -)-e(^W-'0/T- (^9) 

Here h is the chemical potential and B is the external constant field. The 
functions E±(k) are 

E+{k) = E{k)E^{k), E^{k) = exp . (40) 
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We introduce also the function 



Fir,v) :=l + ^7"^(e^'"' + e-'P''), 



where 



J — 1 + e 



2B/T 



(41) 



(42) 



It is worth mentioning that 7 € [l,oo) for any real external field B. At the 
point 7 = 1 (i? = —00) one has F{1; rj) — 27rA(?7) where A(77) is the 27r-periodic 
delta-function. 

Now we are ready to formulate the main results. 

For the correlation function of the canonical field operators of the Hubbard 
model at C/ = 00 we obtain the following representations: 

TT 

— TT 

det (j + 7 Q{t]) + 7?,(-)) - det (j + 7 2(77) 

TT 

(c,,,(i)cS,^(0))(^)=e^*(''-^)^ J dr^Fir,v) 

— TT 

det (i + 7 2(77) -7 7tW(7?)) + (G- 1) det (1 + 72(77) ) , (43) 



(44) 



where X is the unit operator on the interval [— 7r,7r], and 
2(77)^V(77)-l^G7^(-). 



The integral operators V(77), iZ^ \ (77) (the latter two are of rank one) 
possess kernels 



V{v\k,k') = 



£+{r]\k)£-{k') ~e^{k)e+{k') 
27rtan(^) ' 



7^(-)(fc,A;')-:^^-(fc)i?-(fc'), 



2tt 



n(+Hv\k,k') = — ^ -^+(77|A:)£+(77|fc'). 

7r(l — cos 77) 



(45) 



The representations for the correlators {cl |(i)cQ |(0))^"^^ and (c^ j,(Oco |(0))^"^^ 
can be obtained from (Esh by inversing the sign of the external field, B —B. 
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In the equal-time case, t = 0, the integrals over 77 in (^3|) can be taken 
explicitly. For t = and .t 7^ one has 



= det (l - l)v + - det (^1 + (7 - 1)? 
where the kernels of integral operators v and f are 



sm X 



27rtan(^) 



r{k, k') = \A9(fc)e-"^ ^/d{k' 
For i = and a; = one has 

7r 



(46) 



(47) 



(48) 



(49) 



in agreement with the relation (pOD. 

Consider now the temperature correlation functions of the density operators. 
For these correlation functions the following representation is valid 



and therefore, e.g.. 



(T) 



1 

2^ 



dry (7; 77) 



1 — cos 77 



det (i + 7^^(77)) - det (i + jQir])^ 



(50) 



where 



:^i + 2^A(77)^, 



$(Ta)(^. ^) ^ $(i.T)(^. ^) ^ „2_i + 27rA(77) 



7-1 



(51) 
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and A (77) is the 27r-periodic delta- function. The integral operator h({T]) possesses 
the kernel 



U{r]\k,k') = 



27rsin(^^) 



(52) 



and the integral operator Q{t]) is defined in (Q. Let us note that the contribu- 
tions to the correlators (^0|) containing the delta- function A{'r]) in the quantities 
$("'^)(7;77) admit more explicit form. Indeed, only the Fredholm minors up to 
the second order contribute to the correlators in these terms. Let us denote this 
contribution (up to a numerical factor depending on spins) as Q{x, t; h, B). One 
has 



G{x,t-h,B) 



1 — cos 77 



det Z + 7^Y(r/) - det X -I- 7 Q(7/) 



i)=0 



^ite{k)—ixk 



2^ / d/^e 



^ J dfci?(fc), 



,its{k)-- 



(53) 



Other contributions to the correlators (terms independent on 77 in <i>("'^^ (7; 77)) 
could not be simplified in the general case since the integral over 77 cannot be 
taken in the closed form. One can, however, consider some combinations of the 
correlators (|50| ) in which these terms cancel; e.g., 

((n.,T(i)+n.,x(0)(^io,T(0) + rioa(0)))^^^ =^(a;,t;/i,S). (54) 

In the equal-time case the integral over 77 in ( |50| ) can be taken for all the 
contributions. Remarkably that in the case t — and x ^ the correlation 
functions are determined by the function G{x, t\ h, B) only. In this case one has 



— / dfc7?(fc), 

2tt J 



— ixk 



where 



0(TT)^J_^ ^(Ti) ^ ^UT) ^ 2_i^ 



7-1 

7 



(55) 



(56) 



In the case t — and a; = not only the function Q(x,t] h, B) contributes to 
the correlators. The result is 
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— TT 

K,T«^a>*^^ = K,i«-,T>^^^ = 0. (57) 

Let us note that the obtained resuhs for the correlation functions (^ ) and 
( pO| ) has the proper "one-component hmit". Indeed, in the hmit B ^ — oo, 
ft, — > — cxD (h — B = ho is fixed) one has a free fermion model (of fermions with 
spin I). Since 7 = 1 in this limit, the representations ( ^ ) and ( ^ ) become in 
fact trivial reproducing the well-known results for the correlators of free fermions 
on the lattice. Also for the only non-vanishing correlator of density operators, 
(na;_|(/;)no.|(0))^^\ one gets in this limit the expression for the correlator of 
third local spin components in the XXO spin chain 1 23 Uu . 



In conclusion, we would like to stress that the integral operators involved into 



the representations ( [43| ) and (50) are of the form usual for integrable models, 
i.e., they are "integrable integral operators" This fact is important 

for constructing the corresponding matrix Riemann-Hilbert problem and for 
deriving integrable partial differential equations for the correlators. This, in 
turn, will make possible the evaluation of different (e.g., large time and distance) 
asymptotics of the correlators considered. 
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